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Abstract 

We study how the combined observation of dark matter in various direct detection 
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V/^ ■ experiments can be used to determine the phenomenological properties of WIMP dark 

t^ . matter: mass, spin-dependent (SD) and spin- independent (SI) scattering cross section off 

nucleons. A convenient choice of target materials, including nuclei that couple to dark 

^— V \ matter particles through a significantly different ratio of SD vs SI interactions, could break 

Cn ' the degeneracies in the determination of those parameters that a single experiment cannot 

discriminate. In this work we investigate different targets that can be used as scintillating 
bolometers and could provide complementary information to germanium and xenon de- 

^^ , tectors. We observe that AI2O3 and LiF bolometers could allow a good reconstruction of 

^ \ the DM properties over regions of the parameter space with a SD scattering cross section 

as small as 10~^ pb and a SI cross section as small as 5 x 10~^° pb for a 50 GeV WIMP. 
In the case of a CaW04 bolometer the area in which full complementarity is obtained is 
smaller but we show that it can be used to determine the WIMP mass and its SI cross 
section. For each target we study the required exposure and background. 
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1 Introduction 

The detection and identification of the dark matter (DM) is a challenging goal that is currently 
being pursued by a large number of experiments around the world using different techniques. 
DM can be searched for directly (attempting to observe its scattering off nuclei in a detector) , 
indirectly (looking for the products of its annihilation or decay), and in particle colliders 
(which explore the nature of physics at the TeV scale, the typical scale for many models 
of particle DM). Such variety of strategies is sensitive to DM candidates with very different 
properties, allowing the exploration of a wide range of particle models. Among the various 
possibilities, a generic weakly-interacting massive particle (WIMP) is well motivated since it 
predicts quite naturally a value for the DM thermal relic density that is of the same order of 
magnitude as the observed cold DM relic abundance. 

Direct DM detection is currently undergoing an exciting situation, with several experi- 
mental collaborations (using different targets and techniques) reporting potential signals of 
WIMP DM that have, nevertheless, not showed up in other detectors. In particular, an 
annual modulation in the detection rate was observed by the DAMA collaboration using 
Nal(Tl) as target [1] and was later confirmed by the extended experiment DAMA/LIBRA 
(2] reaching a statistical significance of 8.9 cj. A new experiment, ANAIS, projected to be 
carried out at the Canfranc Underground Laboratory with up to 250 kg of Nal(Tl) would test 
this observation in a model independent way OHj. Also, the CoGeNT collaboration (using a 
germanium target) observed an irreducible excess in their data [5] that, if interpreted in terms 
of WIMPs, would correspond to a very light particle, with mass in the range 7—12 GeV, 
and a large elastic scattering cross section, of order 10~^ pb. Furthermore, hints of an annual 
modulation in the CoGeNT experiment were also claimed [6] but with a limited statistical 
significance of 2.8 a in the first year data. Moreover, data presented by the CRESST collab- 
oration (which uses CaW04 as target and beta/gamma background rejection) also displayed 
an excess that could be compatible with light WIMPs • Several analyses have investigated 
the compatibility of these three signals [HJ [9l [TOl [TTl [12] , although this only seems possible if 
extreme assumptions are made for the different uncertainties, such as the inclusion of large 
quenching factors or channeling effects. 

These observations are challenged by the negative results obtained by other experimental 
collaborations. Most notably, CDMS II [ISIIII], XENONIO [E], XENONIOO [El [IT], SIM- 
PLE [E], KIMS [E] and a combination of CDMS and EDELWEISS data [20] have set upper 
bounds on the spin-independent (SI) part of the WIMP-nucleon cross section that are in 
strong tension with the regions of the parameter space compatible with WIMP signals in the 
DAMA/LIBRA or CoGeNT experiments. Several alternatives to ordinary WIMP DM have 
been proposed in the literature, trying to account for this discrepancy. Moreover, CDMS II 
has recently searched for annual modulation in a reanalysis of their data and did not observe 



any indication for it |21j . This further strengthens its incompatibihty with CoGeNT since 
both use the same target, leaving no room for model dependence affecting the comparison. 
Similarly, the spin-dependent (SD) contribution to the WIMP-nucleon cross section is also 
constrained from negative results in direct detection experiments. The leading constraints 
are due to XENONIOO [I71[22l[23] (SD cross section with neutrons, o-"^^'") and COUPP [21], 
PICASSO [25], SIMPLE [18] and KIMS [E] (SD cross section with protons, ct'^^'P). 

In the near future more sensitive experiments are going to continue probing the DM 
parameter space. This will allow us to clarify the current situation regarding light WIMPs 
and also explore DM candidates with smaller interaction cross sections. In particular, some 
of the existing experiments are involved in the upgrading of their detectors, moving towards 
several hundreds kilograms of target material or even reaching the 1 ton scale. This is the 
case of the XENONIT |26] and LUX [27] collaborations, with a xenon target, SuperCDMS 
[28] . using germanium, COUPP, with a C3FI target, and ArDM, using liquid argon [29]. On 
a similar timescale, the EURECA |30j consortium has plans for a 1 ton scale experiment, 
able to operate different types of cryogenic detectors consisting of Ge bolometers measuring 
heat and ionization, and CaW04 scintillating bolometers measuring heat and light, although 
other scintillating targets could also be accommodated. The CRESST [31] and ROSEBUD 
|32] collaborations have worked on R&D and tested other potential scintillating targets, 
focusing especially on targets that contain nuclei with enhanced sensitivity to SD interactions 
and low mass WIMPs. According to the characterization and performances of the different 
scintillating bolometers and their complementarity with Ge and Xe for the determination of 
the DM properties, some of these materials could be incorporated as additional targets for 
EURECA, probably in the second phase. 

If DM is detected, the use of different targets is crucial, as it can serve to unambiguously 
determine some of the WIMP properties (e.g., its mass and interaction cross section off 
protons and neutrons), thus helping to discriminate among the various WIMP candidates. 
This idea was applied to the case of the COUPP experiment in Ref. [33], emphasizing the 
relevant role of targets which are sensitive to SD WIMP-nucleus interactions and showing 
how detection in two complementary targets (in that case C4F10 and C3FI) could allow a 
better measurement of the WIMP couplings. The idea of target complementarity has later 
been applied to the determination of the WIMP mass and cross section from different DM 
experiments [34:\ [35] and the relevance of targets sensitive to the SD cross section has been 
analytically studied in Ref. 



In this paper we focus on the determination of the phenomenological parameters of a 
generic WIMP, namely its mass, m^, the SI contribution to the WIMP-nucleon cross section, 
o""^^, and the SD component, a^^. In doing so we take into account all known sources of 
uncertainties, including those in the astrophysical parameters describing the DM halo, but 



also considering uncertainties in the nuclear form factors along the lines of our previous 
analysis in Ref. [37] . We follow an approach based on Bayesian statistics that allows us to 
incorporate these uncertainties in a consistent way, thereby quantifying their effect. 

With these tools, we carry out a systematic study of the performance of various targets 
used in direct DM searches. We first show that 1 ton experiments based on germanium and 
xenon might be unable to determine some of the WIMP parameters (in particular failing to 
measure the SD component of the cross section). We also show that although C3FI is a good 
alternative, the complementary capability of COUPP is limited by its detection technique 
since it does not provide information about the recoil spectrum. We then turn our atten- 
tion to other possible experiments and study how their use in combination with germanium 
and xenon can serve to unambiguously determine my^, cr , and a in certain scenarios, a 
situation that we define as complementarity. We concentrate on the scintillating bolometers 
currently used by CRESST (CaW04) and on those characterised by the ROSEBUD collabo- 
ration (AI2O3 and LiF), which contain nuclei which are sensitive to the WIMP SD coupling 
and are also optimal for searches for low mass WIMPs. We show how for a certain range in 
the phenomenological parameter space, AI2O3 and LiF can provide a good complementary 
measurement which allows the degeneracy to be disentangled in the SI and SD contributions 
to the total cross section. This is generally the case when the detection rate in germanium 
and xenon is dominated by the SI component. On the other hand, CaW04 is a more con- 
venient target when the rate in germanium or xenon is mostly SD (since tungsten is heavier 
and more sensitive to the SI component) and can be used to provide a good measurement of 

The paper is organised as follows. In Section[2] we summarise the physics of direct DM 
detection and the determination of WIMP parameters from the observed number of events 
and energy spectrum. We describe the reconstruction method used in this work, based on 
Bayesian statistics and show explicitly that conventional targets which are more sensitive 
to the SI WIMP-nucleus cross section (such as Xe and Ge) do not, in general, allow for a 
complete reconstruction of all the WIMP parameters. We also show how the situation changes 
when a target that is sensitive to SD interactions (using a hypothetical 1 ton version of the 
COUPP experiment). In Section[3]we study various possible targets used in the CRESST and 
ROSEBUD experiments and analyse the conditions under which complementarity is achieved. 
The conclusions are left for Section^ 



2 Determination of WIMP properties from direct detection 

We start by summarising the basic expressions that describe the WIMP rate in a direct DM 
detection experiment [S5] (for a recent review see Ref. [SU]). The differential event rate for 



the elastic scattering of a WIMP with mass m^ off a nucleus with mass 7717V and a detector 
mass Mdet is given by 

vf{v) {v,ER)dv, (2.1) 



where po is the local WIMP density, f{v) is the WIMP speed distribution normalized to unity 
and velocities are expressed in the detector reference frame. The integration over the WIMP 
velocity is performed from the minimum WIMP speed needed to induce a recoil of energy En, 
Vmin = \/ {"iTT-N Er) / (2^^^^) where /iAT is the WIMP-nucleus reduced mass, to a escape velocity 
Vesc also in the detector reference frame, above which WIMPs are not gravitationally bound 
to the Milky Way. The WIMP-nucleus differential cross section dawN /dEn is computed 
from the Lagrangian that describes the interaction of a given WIMP with ordinary matter 
and encodes the particle physics input. The total event rate is then calculated by integrating 
the differential event rate over all the possible recoil energies, 

R= I dEn / vf{v) {v,ER)dv. (2.2) 

The window for WIMP searches is selected from a threshold energy, Ex-, to a maximal recoil 
energy, Emax, and depends on the specific experiment. The total number of DM recoils is 
N = Rt, where t is the live time of the experiment. In the following, we will also refer to the 
exposure, defined as e = M^ett- 

In general, the WIMP-nucleon scattering cross section can be separated into a SI and a SD 
contribution. The total WIMP-nucleus cross section is calculated by adding these coherently, 
using nuclear wave functions. The differential cross section thus reads 

-dE^ = 'M4PK'' ^57(^«) + ^o Fsd{Er)), (2.3) 

where <Tq ' and cJq ' are the SI and SD components of the WIMP-nucleus cross sections 
at zero momentum transfer, and the form factors -FJj ^j^IEr) account for the coherence loss 
which leads to a suppression in the event rate for heavy nuclei in the SI and SD contributions. 
See Ref. [10] for a complete description of these prescriptions. 

The WIMP-nucleus cross section at zero momentum transfer can be written as [401 HT] 

SD,N 32fij^Gp r c- , Cl2 Z'^+lA /o /^^ 

o-Q = ^^— ^ [apSp + a„5„,J ( — — 1 , (2.4) 

where Sp and 5„ are the expectation values of the total proton and neutron spin operators; 
/p) fn and fflp, an are the effective WIMP couplings to protons and neutrons in the SI and 
SD case, respectively; Gp is the Fermi coupling constant, and J is the total nuclear spin|^ 



^Notice that for simplicity we have not included here a possible vector coupling (corresponding to non- 
Majorana DM particles), which would lead to an extra contribution in the expression for a^ ' . 



In the following we will assume that fp = fn, so that 



,SI,N _ ( [m\\2^si ^ (2.5) 



where a is the WIMP-nucleon SI cross section and /i„ is the WIMP-nucleon reduced mass. 
In the SD dependent case we a 
WIMP-neutron cross sections as 



In the SD dependent case we can also express CTq ' in terms of the WIMP-proton and 



We will assume a specific relation between the couplings to protons and neutrons, namely 
dp /an = — 1 so that only one parameter a '^ = u '" = a will be needed to describe SD 
interactional (see Section [2rT]) . 

In this paper we follow a phenomenological approach, where no specific particle physics 
model for the DM is assumed. Instead, the WIMP is characterised simply by its mass m^, 
spin-independent and spin-dependent WIMP-nucleon interaction cross sections, a and a . 

If a DM signal is obtained in a direct detection experiment, the observed number of events 
and (if the experiment provides it) the energy dependence of the differential rate, i.e., the 
energy spectrum, can be used to reconstruct the properties of the DM particle [33l [HI BS l [3^ . 
In doing this, it is important to remember that expression ()2.ip is subject to uncertainties 
in the nuclear form factors and in the parameters describing the DM halo. Determining the 
impact of these is crucial to understand the capability of a DM experiment to reconstruct the 
WIMP properties. In particular, astrophysical uncertainties are known to significantly affect 
the reconstruction of the mass and scattering cross section of the DM, see e.g., Refs. |46ll47j . 
Similarly, uncertainties in the SD structure functions can lead to a mis-reconstruction of the 
WIMP mass and SD scattering cross section |37j . 

For a given experimental setup (we use the label a to denote the target) we define an 
energy window for WIMP searches, from a threshold, E^ to a maximum energy -E^ax- We 
divide that energy range into energy bins {Ef,Ef + AE""} with a width AE"". We then 
compute, for a choice of DM parameters, the expected number of events {A?} in each energy 
bin, by integrating Eq. (j2.ip in the corresponding interval for a given live time and adding a 
certain level of background events B^. The specific energy windows, bin size and background 
assumed can be found in AppendixKl for the experimental set-ups considered. 

We consider the quantities {A"} as the experimental information from which we attempt 
to reconstruct the DM parameters. Our analysis is based on the Bayes theorem, which 
determines the posterior probability distribution function (pdf) j»(0|D) of a set of parameters 



^The analysis can be extended to consider a larger parameter space, which would require further study. This 
simplification is motivated by the fact that particle models for DM matter generally predict Ict^^'"! « |a"'^^'^|. 



(for which a prior probabihty is assumed p(©)) from a set of experimental data D, encoded 
in the hkehhood function p(D\@) (or £(©)), 

,_, , p(D|0)p(0) , , 

^(Q|D) = ^ p(Dy ' • (2-7) 

The evidence p(D) in the denominator of Eq. (j2.7p is a function of only the experimental 
data. For our purposes it works as a normalization factor and can therefore be ignored. 
The pdf in Eq. ()2.7p . in principle, depends on the priors p{@) and different choices of priors 
can affect the shape of the final pdf. However, should this happen, it would mean that 
the experimental data are not constraining enough, not being able to dominate the final 
probability distribution. Residual prior dependence can be seen, e.g., in Refs. [Ml SHI [50] . 
The scans over the parameter space that allow us to retrieve the pdf are performed with 
MultiNest 2.9 [511 [52] interfaced with our own code for the computation of the number of 
recoil events and the likelihood. Scans are performed with 20000 live points and a tolerance 
of 0.0001 to reach a good sampling of the profile likelihood (see below) as found in Ref. [49] . 

In our case the experimental data consists of the predicted sets of binned WIMP induced 
nuclear recoils for each target, D = {{Xf})- The parameter space is = {m^, a^^ , a^^) 
and our scans span the following range: m,^ = 1 — 10^ GeV, a^^ = 10^^^ — 10~^ pb, and 
a^^ = 10~*— 1 pb. Logarithmic flat priors are assumed for the three variables. The parameter 
space is further extended with inputs describing the astrophysical and nuclear uncertainties, 
which are considered nuisance parameters. See Appendix[Bl for a detailed description of the 
parameterization of nuclear uncertainties for the different targets. Astrophysical uncertainties 
have been included as in Ref. [37]. Unless otherwise stated, all the scans in this work include 
both astrophysical and nuclear uncertainties. 

The likelihood C{@) is computed for each point in the scan, computing the number of 
recoil events Nf in the i-th. bin for the experiment a, and comparing it with the prediction 
of the benchmark model in the same bin. A" for the same target, assuming that data from 
each experiment follow independent Poissonian distributions, 

^ ^ ^ Aff'(0)^?e-^"(®) 

/:(e)=n^"(e)=nn A-! — • ^^-^^ 

a a i * ' 

Notice that this is equivalent to the product of the likelihoods for each experiment C^. The 
number of recoil events Nf in the i-th bin for each experiment is obtained by integrating 
Eq. (|2.ip between E^ and E'^ + Ai?", for a given live time and including a certain number 
of background events If j l^E"' . 6" indicates the number of background events in the whole 
energy for target a, so that, assuming no energy dependence, 6"/Ai?" is the predicted number 
of background events in each energy bin. We scan over 5" , included in the scan as a nuisance 
parameter. We assume the number of background events 6" follows a Poissonian distribution 



function with a mean value that coincides with the nominal background level expected for 
the different experiments (see Appendix [A|) . 

In the next sections, the results of our scans will be plotted by means of two-dimensional 
plots. When the probability for a subset of the original is considered, one can account for 
the presence of the hidden parameters in two different ways: 

• by marginalizing over them, obtaining the pdf for the j'-th parameter integrating over 
all the others 

p(ej|D) = /"p(0|D) dQi... dOj^i dQj+i dG„; (2.9) 

• by maximizing over them, obtaining the so-called profile likelihood (PL) 

£(9,) = max £{&). (2.10) 

fc)i,...,feij_i,feij+i,fei„ 

The PL is usually more sensitive to small fine-tuned regions with large likelihood, while the 
integration implemented for the pdf allows to account for volume effects. Thus, a parameter 
space characterized by a complicated likelihood function C{@) may result in different pdf and 
PL for the same parameters. In our work, to avoid the inclusion of too many figures, we will 
only represent the results for the PL. In the following sections, we demand closed contours in 
the PL 68% and 99% confidence regions (i.e. a good reconstruction of the DM parameters) 
as a requirement for complementarity. In general, when closed contours are obtained in the 
distribution of the PL for the DM parameters we have also observed closed contours in the 
pdf. 

2.1 Complementarity of direct dark matter experiments 

The detection of WIMP DM in more than one target could provide more information about 
the nature of this particle [331 EH ESI E5]- First, the consistency of the energy spectra 
measured by experiments using different target nuclei would confirm that the events were 
due to WIMP scattering (rather than, for instance, neutron background) [iQ]. Furthermore, 
part of the astrophysical uncertainty can potentially be removed, and in fact the comparison 
of the differential rates could be used to determine the moments of the halo WIMP velocity 
distribution [45^1531 . This is one of the motivations for the multitarget project EURECA. 

The use of targets which are sensitive to both the SI and SD components of the WIMP- 
nucleus cross section might allow us to determine the WIMP couplings to matter in an 
unambiguous way, potentially discriminating between different DM models. In Ref. [33] this 
idea was put to the test by studying two of the possible targets of the COUPP experiment 
(C3FI and C4F10). These two materials are sensitive to different combinations of SI and SD 
WIMP-nucleon couplings and it was shown that this information could be used, in the case 



of a hypothetical detection, to discriminate the hghtest super symmetric particle (LSP) in 
the minimal supersymmetric extension of standard model (MSSM) from the lightest Kaluza- 
Klein particle (LKP) in the universal extra dimension scenario. This is a direction that we 
explore in more detail in this work. In particular, we investigate how well we can reconstruct 
the DM properties assuming an observation in more than one direct detection experiment 
and define complementarity as the capability of a combination of targets to determine these 
properties with a certain precision. 

Let us start by assuming a future detection of WIMP DM in a single direct detection 
experiment. As explained in the former subsection, the phenomenological parameters defining 
the WIMP DM particle can be extracted from the study of the total number of events and 
the differential recoil rate. However, one can easily understand that this procedure does 
not lead to an unique solution and the parameter space is indeed affected by degeneracies, 
since in general the experimental information is not enough to constrain the three unknown 
quantities. 

To illustrate this more clearly, we note that the total number of WIMP recoil events can 
be expressed as 

N = Csia^^ +CsD (2SpVaSD,P + 2SnVaSD,n\\ (2.11) 



where the coefficients Cgj/gD contain the integration in velocities and energies (and a depen- 
dence on the WIMP mass) for a given exposure. 



epo/(v) 



Notice that all the dependence on the astrophysical halo parameters and most of the depen- 
dence on experimental setup, such as target material, energy threshold, energy resolution, 
are contained in them (there is also a dependence on the target material in Sp and Sn). 
As previously stated, we particularise our analysis to the specific case a^^'P = o-^"^." and 
therefore Eq. (|2.1ip contains three unknown quantities: the SI and SD components of the 
WIMP scattering cross section and its mass. Even if we assume that the WIMP mass can be 
determined independently with a reasonable accuracy, we are still left with two parameters 
to reconstruct. Thus, given only one experimental result, the same number of events can 
be accounted for by different combinations of SI and SD couplings. This is illustrated in 
Fig.[Tl where each of the color shaded areas corresponds schematically to the region in the 
{a^^ ,a^^) plane that is compatible with the detection of a certain number of recoils in one 
particular detector. The detection of a WIMP in a second experiment with a different target 
can provide complementary information with which this degeneracy can be partially resolved, 
since changing target implies that also the Cgj/gD coefficients in Eq. (|2.11|) are different. 
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Figure 1: Schematic view of the reconstruction of a^^ and a^^ from the observed number of WIMP 
recoil events in two different DM experiments (orange and green shaded areas, respectively). Even 
assuming a precise measurement of the WIMP mass, a single experiment cannot unambiguously 
determine the SD and SI cross sections. The detection on a second experiment, however, provides 
extra information that may allow to further constrain these two parameters. On the left hand-side we 
display an example in which experiments are complementary, i.e. they intersect leaving only a closed 
region compatible with both data sets. On the right hand-side we illustrate the case in which the two 
experiments are not totally complementary, since the region where they overlap is not bounded from 
below. The nominal WIMP parameters for each case are shown with a dot. 

Fig.[T] depicts two possible situations. On the left hand-side we consider an example in 
which two targets are complementary and allow a good reconstruction of both the SD and SI 
couplings. This is the case, for example, if one target is mostly sensitive to SI interactions and 
the other one to SD ones. It can be seen that the region compatible with both experiments 
is closed, a situation that we call complementary. On the right hand-side we show another 
case in which complementarity is not present since the overlapping region is unbounded 
from below. This is generically the case when the two targets are mostly sensitive to the 
SI coupling, and therefore, is a very common situation. In particular, as we will argue in 
Section[221 this could happen when combining signals from germanium and xenon detectors. 
Notice that this example suggests that in such a case the SI coupling can be relatively well 
reconstructed but only an upper bound can be derived for the SD one. 

The extension of the overlapping region is a result of the various uncertainties. First of 
all, there is an obvious statistical component: if the detected number of events is small, it is 
subject to large statistical fluctuations and therefore the reconstruction of parameters is poor. 
Thus, the larger the exposure of the experiment (or alternatively, the larger the scattering 
cross section), the narrower the reconstructed band is. On top of this, astrophysical and 



10 



nuclear uncertainties further contribute to worsening the reconstruction. This is clearly a 
limiting factor that needs a careful implementation in order to describe a realistic experimen- 
tal situation. Finally, in the examples above we have assumed that the mass of the WIMP 
is known, but in the following the WIMP mass is part of the parameter space we scan over 
and it will be determined along with the rest of the parameters. 

2.2 Combination of signals from germanium and xenon experiments 

Two technologies which have been very successfully applied to direct DM searches are semi- 
conductor cryogenic detectors (e.g., the germanium detectors of CDMS and EDELWEISS) 
and noble liquid detectors (a xenon target in the case of XENONIO and XENONIOO). Both 
techniques show great potential and there are plans to extend target masses up to several 
hundred kilograms or even one ton within the next years (e.g., Super CDMS, EURECA and 
XENONIT). It is therefore conceivable that a future WIMP observation could take place in 
any of these targets. Let us therefore start by contemplating that possibility and assessing 
how well the DM properties can be reconstructed in germanium and xenon targets for dif- 
ferent WIMP scenarios. For concreteness we consider the energy window for WIMP searches 
corresponding to CDMS-II in the case of germanium, XENONIOO in the case of a xenon 
target (see Appendix [A|) . 

We consider the set of WIMP benchmark points from Tabledl which take into account 
various mass choices and include cases in which either the SI or the SD contribution dominates 
the detection rate. We simulate the expected differential spectrum by considering a 1 ton 
experiment with a 30% live time operating for a year, i.e. a total exposure of e = 300 kgyr. 
The expected DM signal is computed as described in the previous section, and implemented in 
the scan as experimental information, in the attempt to reconstruct the DM phenomenological 
parameters. 

Figure[2] shows the PL of the three phenomenological parameters m^, a , and a 
projected onto the three two-dimensional planes {a^^ , m^), {a^^, m^), and {a^^ , o"^^)i for 
a DM particle corresponding to benchmark M-SI. The first row shows the determination of 
parameters after a detection in a germanium detector alone. The yellow dot indicates the 
nominal value of the benchmark point and the circled cross the best-fit point. Consistently 
with the previous discussion, we observe a large degeneracy in the three parameters, with 
a significant uncertainty in both components of the cross section and the mass. As already 
pointed out in a previous work [37] , the reconstructed values of the SD or SI cross section 
show no upper bound for large values of the WIMP mass. Consequently, the 68% confidence 
level contours are not closed, but extend beyond the limits of the plots. Moreover, the scan 
does not manage to reconstruct the correct value of the a . 
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Benchmark Point 


m^ (GeV) 


<y'' (pb) 


a^^ (pb) 


Ngc 


iVxc 


A^CsFI 


M-SI 


100 


10-9 


10-5 


37 (37) 


54 (52) 


42 (32) 


L-SI 


50 


10-9 


10-5 


40 (40) 


64 (62) 


49 (35) 


L-SD 


50 


4 X 10-1° 


6 X 10-^ 


44 (16) 


93 (14) 


560 (16) 


VL-SI 


15 


10-8 


10-5 


29 (29) 


28 (28) 


15(9) 



Table 1: Set of benchmark points used in this work. We consider three regimes for the WIMP 
mass, very hght (VL), hght (L) and medium (M). The label SI or SD in each benchmark indicates 
which component of the scattering cross section dominates the detection rate in germanium and xenon 
targets. For reference we also include the expected number of WIMP recoil events in Ge, Xe and C3FI 
for an exposure e = 300 kgyr, where the number in parenthesis corresponds to the contribution from 
only SI interactions. See Appendix|X] for details on the experimental set-ups. 



The second row in Fig. [2] corresponds to the combination of data from germanium and 
xenon detectors. An improvement is visible since large regions of the parameter space are 
associated with a smaller value of the PL and the best-fit point is now closer to the nominal 
value. However, both the shape and the area of the outer contours remain very similar to 
the case with only germanium. 

Finally, the third row of Fig. [2] considers the combination with data from a hypothetical 
1 ton version of COUPP. Since the target material, C3FI, incorporates fluorine (which has 
a large nuclear spin) it has a larger sensitivity to the SD component. Thus, the allowed 
region is drastically reduced. In particular, the solutions with large SD cross section are not 
compatible anymore with the data, a more stringent upper limit on a is produced (at least 
for the 99% confidence region), which translates into a lower limit on a . Note that the 
experimental setup in COUPP does not allow the recoil energy to be measured and therefore 
does not provide more information on the WIMP mass. 

Blue contours in Fig. [2] indicate the 68% and 99% confidence regions for the PL of scans 
performed without including the nuclear uncertainties associated with the SD structure fac- 
tors and fixing their parameters to the central value of the range considered in our scans (see 
AppendixlBj) . We observe that the effect of these uncertainties can greatly affect the recon- 
struction of parameters, significantly enlarging the areas compatible with the observation in 
two experiments (see the second row of Fig.[2|). 

These results depend on the benchmark point, but the general conclusion on how different 
experiments combine remains valid. When the number of events increases, either because 
the scattering cross section is larger or because the WIMP mass is closer to the optimal 
sensitivity for these targets, the determination of the DM parameters is better. For example, 
Fig. [3] represents the hypothetical detection of a WIMP with properties determined by case 
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planes for the benchmark point M-SI. The plots in the first row are obtained assuming a detection 
only in a Ge experiment, the second row shows the combination between Ge and Xe targets, and the 
third row corresponds to Ge-l-Xe-t-CaFI. From the inside out, contours are 68% and 99% confidence 
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L-SI in Tabledl in which m^ = 50 GeV. Since the number of recoils is higher, the 68% and 



70 contours in the PL are smaller. In particular, using only data from Ge and Xe detectors 
there is a good reconstruction of the WIMP mass, as well as a significant improvement 
in the determination of a^^ . Notice that even then, the SD cross section is not properly 
determined, as only an upper bound is obtained. As before, after the inclusion of C3FI data, 
the SD cross section is better reconstructed (i.e. lower and upper limits are obtained) at the 



14 



J3 
Q. 



O 




° -5h 



L-SD ' ',^^;?^ 


■°= ..^.-^^ ^- 


i 


i^Fr^ 


< 

■ < 


i 


< 




< 

< 


■! 



S--3 
-9- 
M -4 

° -5 



H 




4 


r . 


^ 


^^^^H ' 




- 


i 


■ 


> 


- L-SD 




^B 


- 


Ge 




^^H 


> , 



-11 



-10 



log(m /GeV) 



log(m /GeV) 



-9 
log(a^'/pb) 




log(m /GeV) 



log(m /GeV) 



-9 
log(aS'/pb) 



-lOl- 



-11 



L-SD 


Ge+Xe+CFjl 


^B 


\ 


II 


\ , " 




I -3 
CO _4 

° -5 



■ 




^^^- 


- 


i 


■1 


W 


- 


■ 




■ L- 


SD 






Ge 


tXe+CFjl 







-11 



-10 



log(m /GeV) 



log(m /GeV) 



log(aS'/pb) 



Figure 4: The same as in Fig.[2l but for tire benclrmark point L-SD. 



Yo confidence level, with the upper limit being valid also at the 99% confidence level. A 
good reconstruction of a has also the consequence of providing a lower bound for a^^ and 
reducing the uncertainty on m-^^. Both quantities are now determined with an uncertainty of 
approximately one order of magnitude, an extremely interesting scenario made possible by 
the nice interplay between the three experiments. 

The importance of nuclear uncertainties is evidenced when we consider the combination 
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of the signal in Ge and Xe. For a fixed SD form factor it is not possible to provide a good 
fit to both signals with a SD-dominated energy spectrum. This produces an upper limit on 
a'^^. On the other hand, when also the signal in CF3I is included, selecting only points with 
a Sl-dominated rate, the details of the SD form factor become irrelevant again and we see no 
difference between the blue and the black contours. 

The last example considered is given in Fig.U corresponding to the case L-SD, in which 
the SD contribution to the number of DM events dominates. As in the previous cases, the 
combination of data from Ge and Xe is not enough to determine the SI and SD contributions 
(only a combination of both). Including CF3I data has the effect of excluding small values of 
a (contrary to what happened for M-SI and L-SI). Moreover, small SI cross sections also 
produce a slightly worse fit to the data now. Thus complementarity could be obtained for 
the 68% confidence region but only if we do not include nuclear uncertainties. 

To sum up, although 1 ton experiments based on germanium or xenon have an excellent 
discovery potential for WIMP DM, they are mostly sensitive to the SI part of the WIMP- 
nucleus interaction. As a consequence, the combination of data from both detectors would 
generally not allow the determination of the DM parameters in an unambiguous way. If we 
look for complementary information, the obvious alternative is a target that is sensitive to the 
SD component. CF3I is certainly a good choice in this respect, since it incorporates fluorine 
^^F, which has an unpaired proton and a large nuclear nuclear spin J. However, COUPP 
would not provide information about the recoil spectrum, which limits its complementarity 
capability. Even considering a 1 ton phase the degeneracies in the reconstruction of DM 
parameters would not be completely removed. 

3 Complementarity of bolometric targets 

A much more appealing situation could come from an experiment which is sensitive to the SD 
WIMP-nucleon interaction and in which the energy spectrum can be accurately determined, 
as well. This can be the case of certain scintillating bolometric targets that have been 
developed and studied by the ROSEBUD [5l] and CRESST collaborations and are planned 
or could be used by the EURECA project. We will consider the following targets: 

• CaW04: This is a material which has been tested and used as a WIMP target by ROSE- 
BUD [25j and it is the current target in CRESST [7]. The sensitivity to SD interactions 
comes from tungsten, for which the isotope ^^^W (with a natural isotopic abundance 
of 14.3%) has an unpaired neutron. However, tungsten is a very heavy material and 
therefore it is mostly sensitive to the SI component of the WIMP interactions. For 
this reason, we should expect compatibility with experiments based on Ge or Xe to be 
limited. Notice that our analysis does not include the possibility (shown by CRESST) 
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Benchmark Point 


m^ (GeV) 


<y'' (Pb) 


a^^ (pb) 


^CaW04 


A^AbOs 


iVLiF 


M-SI 


100 


10-9 


10-5 


27 (27) 


14(4) 


25(2) 


L-SI 


50 


10-9 


10-5 


30 (30) 


17(5) 


35(2) 


L-SD 


50 


6 X 10-1° 


4 X 10-4 


20 (18) 


505 (3) 


1295 (1) 


VL-SI 


15 


10-8 


10-5 


8(8) 


23 (18) 


29 (11) 



Table 2: Total recoil events for the set of benchmark points expected on each of the bolomctric 
targets considered in this work for an exposure e = 300 kgyr. 



of disentangling O, Ca and W nuclear recoils. 

AI2O3: Sapphire is a very promising DM target because it is sensitive to low mass 
WIMPs (Al and O are both light nuclei). It is also sensitive to SD interactions (^"^Al 
has 100% isotopic abundance and J = 5/2) and recent tests indicate very high light 
yields |561 157] and a particle discrimination threshold below 10 keV seems to be possible 

ESI. 



• LiF: This target is also sensitive to low mass and SD interactions (Li and F are light 
nuclei, ^^F has 100% isotopic abundance and J = 1/2, ^Li has 92.5% isotopic abundance 
and J = 3/2). However, up to now, low particle discrimination thresholds have not 
been obtained and more R&D is needed on this target to be used for DM searches. 

In this section we determine the conditions under which the different scintillating bolo- 
metric targets can provide good complementary information when combined with germanium 
and xenon detectors. As in the previous section, for each benchmark and for each detector 
target, we simulate the number of recoils predicted in the different energy bins. This consti- 
tutes our set of observables. Then, we perform a scan on the three-dimensional parameter 
space {m^^ o"^^, a^^) and apply the Bayesian inference method to determine the reconstruc- 
tion of these quantities. Astrophysical and nuclear uncertainties are included as described 
above. 

We start by considering the same exposure for all the experiments, e = 300 kg yr, and 
study the same benchmark scenarios as in the previous section. In Table[2] we indicate the 
number of recoil events for each of the bolometric targets. Initially we assume zero background 
for the three targets and postpone the discussion about the influence of the background level 
to the end of this section. In Figs.O [6l [7] and [8] we represent the PL for reconstructed 
DM parameters in benchmark points M-SI, L-SI, L-SD, and VL-SI, respectively. The black 
contours correspond to the combination of Ge and Xe data with a bolometric target: CaW04 
for the upper row, AI2O3 in the middle and LiF for the lower row. For comparison, the blue 
contours illustrate the results when only germanium and xenon detectors are used. 
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Figure 5: Profile likelihood for the DM parameters in the {a , n^x); ('^" ' '"-x) 
planes for the benchmark point M-SI after the combination of data from a Ge detector, a Xe detector, 
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is e = 300 kg yr for all the experiments. From the inside out, contours enclose 68% and 99% of 
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Figure 6: The same as in Fig. [5] but for the case of L-SI. 
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Figure 7: The same as in Fig. [5] but for the case of L-SD. 
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We observe how in general AI2O3 and LiF provide good complementarity with germanium 
and xenon in points for which the detection rate in the latter is dominated by SI contributions. 
This is the case of benchmark points M-SI and L-SI. The area in the parameter space that is 
compatible with the simultaneous observation of a WIMP in the three experiments becomes 
much narrower than the case with only Ge and Xe. Although the WIMP mass cannot be 
completely reconstructed in some of these cases (e.g., in M-SI), a lower constraint is generally 
obtained for a (which would suffice, e.g., to discriminate the observation from the case 
a = of scalar DM). In the examples for which the number of DM events is larger, such 
as in L-SI, the extra information on the recoil spectrum results in a better determination of 
the WIMP mass and consequently leads to closed contours in the three quantities. This is 
the situation that we describe as complementary. Although in these cases data from CaW04 
does not allow a complete complementarity, an upper bound can be obtained in the value of 
a^^ which, in turn, leads to a lower bound on a^^ and a good determination of the latter 
(see, e.g., the case of L-SI in the top row of Fig. [6]). 

On the contrary, CaW04 performs better for benchmark points where SD contributions 
dominate the rate for germanium and xenon. Since tungsten is a heavier nucleus than both 
germanium and xenon, it is more sensitive to a . This is, e.g., what happens in L-SD, 
where the reconstructed contours are closed and the reduction in the best fit areas is very 
significant. In this kind of points, AI2O3 and LiF do not perform that well since they are 
mostly sensitive to the SD component. In any case, they can be used to obtain a lower bound 
on a^^ which clearly rules out the possibility a^^ = 0. 

It should also be noticed that in the case of the benchmark point M-SI represented in 
Fig. [5] the reconstruction of the DM mass has no upper bound. This is a generic feature for 
heavy DM particles, due to the fact that the spectrum becomes flatter and the fit to the mass 
is more inaccurate. The combination with bolometric targets does not improve this situation 
significantly. For this reason, we expect a worse reconstruction for DM particles heavier than 
in M-SI. 

Very light WIMPs, on the other hand, constitute an interesting possibility that can also 
be explored with the aid of bolometric targets. In particular, as we said above, both AI2O3 
and LiF can be sensitive to low-mass DM particles. In Fig. [8] we show the case of a WIMP 
with m-f, = 15 GeV and scattering cross section as in VL-SI. Data from germanium and xenon 
are enough to determine the WIMP mass rather accurately, but a large uncertainty remains 
in both a and a as can be observed in the blue contours. Both AI2O3 and LiF remove 
significantly this degeneracy and in the case of LiF we even obtain an inner contour around 
the correct value of a^^ . In the case of CaW04 only the SI component can be determined 
but no further information on a^^ is obtained. 
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Figure 9: The dark blue region represents tire area of the (o"' , o"' ) plane for which complementarity 
is obtained for CaW04, AI2O3 and LiF in the case of a WIMP mass of m^ = 50 GeV. The light blue 
region is a measure of the uncertainty of our grid scan (see text). The black line corresponds to 
the upper constraint obtained from XENONIOO data (using the BonnA result for the SD structure 
function [SH]). Points above the dashed line predict more than 4.9 events in XENONIOO and are 
therefore excluded by the recent experimental result [T7] following a Feldmans-Cousin method [23] . 
For reference, points along the solid lines predict 2 events in XENONIOO. 

So far we have observed that some targets perform better than others (in terms of com- 
plementarity) in certain regions of the parameter space. In particular, for each of the targets 
we can determine the regions in the parameter space for which complementarity is obtained. 
In Fig.[9]we indicate the area (dark blue) of the (o"^^, a^^) plane for which we obtain closed 
contours in the reconstructed PL of the three DM parameters (m^, a"^^, a^^) for each of the 
bolometric targets and a WIMP mass of m^ = 50 GeV. In order to obtain this region we have 
performed a grid scan in the SD-SI plane for which the separation among points is limited 
by computing time. The light blue region separates points leading to complementarity from 
the nearest which do not, hence indicating the resolution of our grid scan. The areas are 
different for the three bolometers studied. Consistently with the individual examples that 
were analysed previously, we observe that the areas for AI2O3 and LiF are larger as com- 
pared with the area for CaW04 , and shifted towards smaller values of a . For CaW04 total 
complementarity only occurs for a small region, very close to the current upper constraint 
by XENONIOO. However, as we showed in various examples, it helps in the determination 
of the WIMP mass and a^^ . Needless to say, a reduction in the exposure would result in a 
shift of the complementarity areas towards larger values of both a and a . 

For completeness, we have determined the highest level of background for which comple- 
mentarity is attained as a function of the exposure for each bolometer in a given benchmark 
point, and represented the results in Fig.llOi We see how, for zero background, the exposure 
for the three bolometers can be reduced. For AI2O3 and LiF this reduction can be very 
significant. For example, in the benchmarks that we have chosen in Fig.llO^ it suffices to 
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Figure 10: The solid line represents the maximum level of background for which complementarity is 
attained as a function of the exposure for each of the bolometric targets for representative benchmark 
points ((cr^^cr^-^) = (6 X 10-i",4 x IQ-**) pb for CaW04, (10-^,5 x IQ-^) pb for AI2O3, and 
(10"^ 10^5) pb for LiF, from left to right) with m^ = 50 GeV. The expected total number of WIMP 
recoil events is indicated in the top horizontal axis. Dotted lines denote the number of background 
events per bin. 

have a clear signal for DM in AI2O3 or LiF, even with a very reduced number of events (to 
determine the lower bound we considered that 1 event in a background free experiment is 
statistically significant). From these results we can also conclude that the complementary 
areas of Fig. [9] would not shrink significantly if a moderate background is included. We have 
to bear in mind that this computation was carried out assuming a flat background, for which 
a DM signal can be easily distinguished. This assumption can be considered equivalent to 
the estimate of the sensitivity of an experiment that does not surpass this background level 
in any bin of the energy window, independently of the background dependence on energy. 



4 Conclusions 

In this work we have investigated the determination of DM parameters from a combined 
use of different targets in direct detection experiments. More specifically, we study the 
determination of the WIMP mass and the SI and SD components of its scattering cross section 
off nucleons. We apply the method to a set of bolometric targets developed and characterised 
by the ROSEBUD collaboration (AI2O3 and LiF) or in use in CRESST (CaW04). We 
investigate the conditions under which the DM parameters can be obtained unambiguously 
when combining data from them and germanium and xenon experiments. In doing so, we take 
into account astrophysical uncertainties in the DM halo properties and nuclear uncertainties 
in the SD structure functions. 

We first show how one ton scale germanium and xenon targets, which might excel in 
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providing the first measurement of a WIMP, miglit not be able to measure all the DM 
properties. This is due to a degeneracy in the SI-SD plane that might be unresolved since both 
targets are mostly sensitive to the SI component of the WIMP cross section. Interestingly, SD- 
sensitive targets might provide extra information with which this degeneracy can be broken. 
We have studied the combination of germanium and xenon based targets with the one used 
by the COUPP collaboration (CF3I). The presence of fluorine improves the reconstruction of 
DM parameters, however, COUPP does not provide information about the recoil spectrum, 
which limits its complementarity capability. 

The situation is much more interesting for experiments which are sensitive to the SD 
WIMP-nucleon interaction and which can provide an accurate measurement of the energy 
spectrum, such as the scintillating bolometric targets developed and studied by the ROSE- 
BUD and CRESST collaborations. We observe how the combined use of these detectors with 
germanium and xenon experiments can break the degeneracies in the determination of DM 
parameters and provide a good reconstruction of the WIMP mass and SI and SD scattering 
cross section. In particular, the inclusion of CaW04 can lead to a better determination of the 
WIMP mass and SI cross section. Nevertheless, since its interaction with DM is dominated 
by the SI contribution, it only leads to a complementary results in a small window of the 
parameter space in which the rate in germanium and xenon is mainly due to SD interactions. 
This area is in fact very close to the region already excluded by XENONIOO. On the other 
hand, AI2O3 and LiF (being more sensitive to the SD contribution) can be complementary 
targets to germanium and xenon in regions of the parameter space where the rate in the 
latter is dominated by SI contributions. This can happen for values of the cross section as 
small as a^^ ^ 2 x 10~^° pb and a^^ > 10~^ pb for a WIMP with a mass m^ = 50 GeV. 

In some regions of the DM parameter space the exposure can be reduced to approximately 
50 kg yr for AI2O3 and LiF without loosing complementarity under the assumption of zero 
background. Finally, we investigate the effect of the background and observe that comple- 
mentarity can be achieved for large exposures with a background level as large as 10~^ kg~^ 
day~^ keV~^ for AI2O3 or 10~^ kg~^ day~^ keV^^ for LiF. In the case of CaW04 a larger 
exposure and smaller background level is required. 

A Experimental features of the different targets 

In this Appendix we detail the experimental parameters that has been used for each detector. 
We have considered natural abundances of the different isotopes for a given target. 

We divide the energy window sensitive to recoils (from Ej- to Emax) into Ni^n evenly 
spaced bins with a size of IS.E = 5 keV, with Niyi^ = {E^^x ~ E^)/ ^E. The estimated 
number of events in each bin is kept as a decimal number, without rounding it to an integer. 
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Target 


{Et, Emax} 

(keV ) 


aiE) 
(keV) 


AE 
(keV) 


Background 
(kg-May-^keV-i) 


Ge 


{10, 100} 


^(0.3)2 + (0.06)2 E/keY 


5 


4 X 10-s 


Xe 


{8.4, 44.6} 


0.6^E/keY 


3.64 


4 X 10-9 


C3FI 


{10, 200} 


- 


- 


4.1 X 10-s 


CaW04 


{10, 100} 


5% FWHM 


5 




AI2O3 


{10, 100} 


5% FWHM 


5 




LiF 


{10, 100} 


5% FWHM 


5 





Table 3: Energy range, parameterization of the resolution and background for the different 
detectors considered in this work. The background level for the first three targets is inspired 
on estimates for SuperCDMS [60], XENONIT [6ll[62], and COUPP [63], that we consider 
energy-independent for simplicity. 



This is, in principle, not physical and might overestimate the ability to discriminate different 
spectra. However it allows us to neglect the dependence of the reconstructed parameters 
on the particular realization chosen for the nominal number of events in the different bins 
(see Ref. [M])- This is an important source of uncertainty that has to be taken into account 
when dealing with real data, but we decide to neglect it here in order to study and estimate 
complementarity in a scenario uncompromised by statistical fluctuations, and also because, 
as found in Ref. [63] the relevance of those statistical fluctuations decreases as the number 
of experiments increases, thus we assume that coverage is good enough when dealing with 
signals from three detectors. Notice that, in the case of LiF, the thresholds obtained to this 
date are far from the 10 keV value used in this work. Additional R&D is needed on this 
target before using it for DM detectors. 

In our analysis the energy resolution of the detector a is included as a convolution of the 
differential rate with a Gaussian function with a variance 172 which depends on the recoil 
energy and on the particular experiment (see Table[3|). 



B Uncertainties in the spin-dependent structure function 



We have incorporated nuclear uncertainties in the SD structure functions for each nucleus. 
Following the procedure described in Ref. [37] we use a three-parameter function to describe 
the various components of the SD structure functions for the various target nuclei. 



5,j(u) = iV((l-/3)e-"" + /3) 



(B.14) 
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Isotope 


N 


a 


/3 


^Li 


0.129-0.077 


2.25 





1^0 


0.133-0.010 


3.16 





19p 


0.144-0.110 


3.03 


0.009 - 


27A1 


0.202 - 0.146 


1.54 





73 Ge 


0.206-0.117 


6.00 - 5.04 


0.04-0.02 


i29Xe 


0.052 - 0.029 


4.66 - 4.20 


0.007-0.001 


i3iXe 


0.025 - 0.017 


5.00 - 4.28 


0.061 - 0.042 


183^ 


0.015 - 0.0005 


3.61 






Table 4: Ranges of the inputs N, a and /? used in the parameterization of the Sn term of 
the SD structure function of the different target nuclei used in the text. A flat probability 
distribution is assumed for each parameter. 

The parameters A^, a and /? are varied in such a way that the family of curves provides an 
envelope for the existing theoretical calculations. These quantities are computed using nuclear 
physics models, and the results may differ, depending on the methodology and the potential 
used to describe the nuclear interaction. For the targets under consideration we extracted the 
predictions from Refs. [65^166] and references therein. Since we have concentrated on the case 
o-p/an = — 1, only the ^n component is important in our analysiso We summarise in Tabled] 
the resulting ranges for N, a and /3 used in our calculations. For some isotopes shell model 
computations of the form factors are not available. In these cases we use the approximation 



Sij{u) = Sij{0)e'' 



^/?2/4 



(B.15) 



which works well in the low momentum transfer regime, and where R is an effective radius 
which is a function of the atomic number A. In these cases a fixed value of a is used and 
/3 = and the only quantity that is varied is A^, which parameterizes the zero momentum 
value of the SD structure function. 
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